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Abstract 

In  the  last  few  years,  a  renewed  interest  in  the  water  gas  shift  (WGS)  reaction  at  low  temperature  has  arisen  due  to  its  application  to  fuel 
cells. 

In  this  work,  a  simulation  of  a  fixed  bed  reactor  for  this  reaction,  which  forms  part  of  a  hydrogen  production-purification  train  for  a  10  kW 
PEM  fuel  cell  using  ethanol  as  the  raw  material,  was  carried  out.  A  commercial  Cu/Zn/Ba/Al203  catalyst  was  employed  and  a  one-dimensional 
heterogeneous  model  was  applied  for  the  simulation.  The  catalyst  deactivation  due  to  thermal  factors  (sintering)  was  taken  into  account  in  the 
model.  Isothermal  and  adiabatic  regimes  were  analyzed  as  well. 

Results  of  the  simulation  indicate  that  the  pellet  can  be  considered  isothermal  but  temperature  gradients  in  the  film  cannot  be  disregarded.  On 
the  other  hand,  concentration  gradients  in  the  film  can  be  ignored  but  CO  profiles  are  established  inside  the  pellet.  Adiabatic  operation  can  be 
recommended  because  of  its  simplicity  of  operation  and  construction.  The  reactor  volume  is  strongly  sensitive  to  the  CO  outlet  concentration 
at  CO  levels  lower  than  6000  ppm.  For  a  lOkW  PEM  fuel  cell,  using  adequate  pellet  size  and  taking  into  account  the  catalyst  deactivation,  a 
reactor  volume  of  0.64 1  would  be  enough  to  obtain  an  outlet  CO  concentration  of  about  7 160  ppm.  This  concentration  value  can  be  handled 
by  the  next  purification  stage,  COPROX. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  clean  and  non-contaminating  characteristics  of  EL  as 
a  fuel  will  depend  on  the  process,  the  raw  material  and  the 
source  of  energy  employed  for  its  production.  If  it  is  obtained 
from  fossil  hydrocarbons,  carbon  dioxide  is  released  to  the 
atmosphere,  regardless  of  the  technology  used;  if  biomass  is 
used  instead  of  hydrocarbons,  carbon  dioxide  is  also  released, 
but  in  this  case  the  plant  consumes  the  carbon  dioxide  pro¬ 
duced  during  its  growth.  Thus,  the  qualification  of  “clean” 
fuel  is  only  true  if  the  raw  material  is  biomass.  Ethanol  is  a 
renewable,  non-toxic  and  easy  to  manipulate  source,  and  has 
therefore  excellent  chances  to  replace  fossil  fuels  [1]. 
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The  new  application  of  EL  as  a  raw  material  for  fuel 
cells  for  mobile  power  sources  (PEM  fuel  cells)  requires 
that  the  anode  inlet  gas  have  a  CO  concentration  lower  than 
10-20  ppm.  Otherwise,  the  anode  is  poisoned  and  the  cell 
efficiency  abruptly  drops. 

Hence,  if  the  H2  is  produced  from  hydrocarbons  or  alco¬ 
hol  reforming,  purification  is  required  in  order  to  reduce  the 
CO  levels  to  cell  requirements.  So  far,  the  most  technologi¬ 
cally  feasible  purification  train  consists  of  a  water  gas  shift 
converter  (WGS)  and  a  latter  step  of  remaining  CO  elimi¬ 
nation.  According  to  Zalc  and  Loffler  [2],  the  WGS  reactor 
is  expected  to  have  the  largest  volume,  and  is  limited  by  the 
catalyst  intrinsic  activity. 

The  reaction 

C0  +  H20  =  C02  +  H2, 

AH°98  =  —4.11  x  1011  erg  mol-1(— 41.1  kJ  mol-1) 
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Nomenclature 

a  activity  (-) 

a0 o  residual  activity  (-) 

A  reactor  cross-section  area  (cm2) 

Bim  mass  Biot  number  (-) 

C  concentration  (mol  cm-3  s- 1 ) 

Cpm  heat  capacity  of  gas  mixture  (erg  mol-1  K-1) 
Ve  reaction  diameter  (cm) 

VQf  effective  diffusion  coefficient  (cm2  s-1) 

F  mole  flow  (mol  s  - 1 ) 

G  mass  flux  (g  cm-2  s  - 1 ) 

AH  reaction  molar  enthalpy  (erg  mol  - 1 ) 

k  reaction  rate  constant  (mol  g  - 1  s  - 1  atm- 1  225 ) 

kd  deactivation  rate  constant  (s-1) 

Ke q  equilibrium  constant  (-) 

km  mass  transfer  coefficient  (cm  s  - 1 ) 

L  reactor  length  (cm) 

n  deactivation  order  (-) 

P  pressure  (atm) 

Pep  Peclet  number  (-) 

A P  pressure  drop  (atm) 

Rq  gas  constant  (cm3  atm  mol-1  K-1) 

Rp  pellet  radius  (cm) 

R\  reaction  rate  (mol  cm-3  s  - 1 ) 

t  time  (s) 

T  temperature  (K) 

uz  Gas  velocity  (cm  s - 1 ) 

v  trial  value  (-) 

V  reactor  volume  (cm3) 

v  CO  conversion  (-) 

y  mole  fraction  (-) 

z  dimensionless  axial  variable  (-) 

Z  axial  variable  (cm) 

Greek  letters 

a  stoichiometric  coefficient  (-) 

p  equilibrium  approach  factor  (-) 

X  auxiliary  variable  N°  1  (-) 

£  error  function  (-) 

£b  bed  porosity  (-) 

0  Thiele  modulus  (-) 

0  effectiveness  factor  (-) 

k  thermal  conductivity  coefficient 

(erg  cm-1  K-1) 

/x  gas  viscosity  (g  cm-1  s-1) 

v  reaction  order  (-) 

0  dimensionless  temperature 

p  gas  density  (gem-3) 

r  dimensionless  reactor  volume  (-) 

co  dimensionless  radial  variable  (pellet)  (-) 

§  dimensionless  concentration  (-) 

\//  auxiliary  variable  N°2  (-) 

£2  radial  variable  (pellet)  (cm) 


Subscripts  and  superscripts 
g  Value  at  bulk  gas 

i  Species  i 

k  kth  iteration 

s  Value  at  surface 

T  Total 

0  Inlet/reference  conditions 


has  been  employed  for  40  years  in  the  industrial  process 
for  H2  production  from  liquid  and  gaseous  hydrocarbons. 
The  role  of  the  WGS  reaction  is  to  increase  the  H2  yield 
and  decrease  the  CO  concentration,  which  is  a  poison  for 
some  catalysts  used  downstream  like  in  ammonia  synthesis 
or  oil  dehydrogenation.  The  WGS  reaction  is  traditionally 
carried  out  in  two  fixed  bed  adiabatic  reactors,  connected  in 
series  with  a  cooler  between  them.  The  first  reactor  operates 
at  temperatures  ranging  from  300  to  500  °C  and  employs  a 
Fe/Cr  catalyst.  The  second  reactor  operates  at  lower  tempera¬ 
tures  (180-300  °C)  in  order  to  displace  the  equilibrium,  since 
the  WGS  reaction  is  exothermic,  and  uses  a  Cu/Zn/Al  cata¬ 
lyst.  Several  studies  that  intended  to  elucidate  the  kinetics 
and  reaction  mechanism  [3-21],  as  well  as  the  deactivation 
by  poisoning  [22],  have  been  carried  out,  Chocron  et  al.  [23] 
have  analyzed  the  mass  transfer  in  a  Cu/Zn/Al  catalyst  subject 
to  deactivation  and  Gonzalez- Velasco  et  al.  [24],  considering 
catalyst  deactivation,  have  proposed  optimal  practices  for  an 
industrial  reactor  operation. 

Recently  a  renewed  interest  in  the  study  of  this  reaction 
has  arisen  due  to  its  application  to  fuel  cells.  Several  works  on 
different  catalysts  more  active  than  the  traditional  one,  which 
can  be  used  for  H2  production  for  a  fuel  cell,  have  been  pub¬ 
lished.  Some  authors  propose  the  doping  of  Cu/Zn/Al  based 
catalysts  [25].  Some  of  them  suggest  the  use  of  the  commer¬ 
cial  Cu/Zn/Al  catalyst  by  adding  O2  to  the  feed  to  oxidize 
the  CO  [26].  Other  authors  have  used  noble  metals  [27,28]. 
Finally,  the  majority  reject  traditional  catalysts  and  propose 
the  use  of  supported  noble  metals,  ceria  being  the  most  com¬ 
monly  used  support  [29-34].  Zalc  et  al.  [33],  after  working 
with  a  1 .5%  Pt/Ce02  catalyst  and  making  an  accurate  analysis 
of  works  in  which  platinum  and  ceria  are  used,  conclude  that 
efforts  should  be  focused  on  finding  commercially  viable  cat¬ 
alysts,  using  lower  cost  metals,  than  those  traditionally  used, 
which  is  a  good  starting  point. 

In  agreement  with  Zalc’s  arguments  [2,33],  a  simulation 
of  a  fixed  bed  reactor  for  the  WGS  reaction,  which  forms  part 
of  a  purification  train  for  a  10  kW  PEM  fuel  cell,  was  carried 
out  in  this  work.  A  commercial  Cu/Zn/Al203  catalyst  doped 
with  Ba  was  used  because  it  showed  a  higher  activity  than  the 
traditional  one.  The  one-dimensional  heterogeneous  model 
was  applied  in  this  simulation,  and  a  parametric  sensitivity 
analysis  was  carried  out  for  some  of  the  process  variables, 
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with  the  purpose  of  finding  criteria  to  minimize  the  reactor 
volume. 


2.  Methodology 

2.7.  Description  of  the  model 

The  heterogeneous  one-dimensional  model  used  in  this 
work  takes  into  account  the  existence  of  concentration  and 
temperature  gradients  inside  the  catalytic  pellet  as  well  as  in 
the  external  film.  Regarding  the  fluid  dynamics,  the  plug  flow 
model  can  be  applicable  if  in  Eq.  (1)  is  fulfilled  [35] 


L  40  vco 
>  - In 
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R 


PeT 


l  —  x 


(1) 


P  A 

A  power  law  intrinsic  kinetics  with  a  term  that  takes  into 
account  the  reversibility  of  this  reaction  is  used  [18] 


where 

Pco2  Pu2  1 


Pco  Ph2o  Kq q 

The  catalyst  deactivation  was  assumed  to  be  only  due  to  ther¬ 
mal  effects  (sintering) 

da 


d  t 


=  -kd(a  -  aD o) 


(4) 


where  a0 0  is  the  residual  activity. The  activity  is  defined  as 
R\(t) 


a(t)  = 


(5) 


*v(0) 

Vi,  k,  n ,  a0 o  and  kd  values  in  Eqs.  (2)  and  (5)  were  experi¬ 
mentally  obtained  using  a  fixed  bed  lab  scale  reactor.  These 
values  are  not  reported  because  of  confidentiality. 

The  system  is  divided  into  two  subsystems:  the  pellet  and 
the  reactor  (or  bulk).  The  pellet  must  be  solved  first.  This  sub¬ 
system  is  a  second  order  non-linear  boundary  value  problem, 
whereas  the  reactor  subsystem  is  a  first  order. 

2.2.  Pellet 

The  Thiele  modulus  is  [36] 


(6) 


In  this  definition  the  surface  and  bulk  concentration  are 
assumed  to  be  equal  since  the  difference  between  them  was 
negligible  (see  Fig.  3). 

The  mass  balance  for  reference  species  CO  is 
—aRy  2  d C 


d  2C 
~dSfi 


=  < 


—OiRy 

3Pef 


Q  d£2 


Q  >  0 


(7) 


In  Eq.  (7)  the  origin  of  coordinates  corresponds  to  the  cen¬ 
tre  of  the  pellet.  It  should  be  noted  that  the  function  including 
the  term  with  the  first  derivative  has  a  singularity  at  Q  =  0, 
so  it  cannot  be  used  at  this  position.  Instead,  the  limit  of  the 
expression  must  be  used. 

Since  the  pellet  can  be  considered  isothermal  (as  it  will 
be  proven  later),  the  transport  properties  are  assumed  to  be 
constant.  Thus,  it  is  possible  to  express  the  species  i  and  T  as 
linear  functions  of  the  reference  species. 


,s  oil  Vef 
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where  the  reference  species  has  no  subscript.The  boundary 
conditions  are 


dC 
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(2) 

d^2 

dC 
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d^2 
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S  =  ^L(c§  -  Cs) 
72ef 


(11) 


Defining  the  mass  Biot  number  (Eq.  (12)),  Eqs.  (10)  and  (11) 
the  dimensionless  boundary  conditions  are  given  by  Eqs.  (13) 
and  (14) 
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Finally,  the  effectiveness  factor  is 
Ry  3 


'Rr 


T]  = 


Rf  4 txRI 


Ry  4nf2~  d Q 


(15) 


.y-rytivp  JQ 


Eq.  (15)  can  be  rearranged  introducing  Eqs.  (6),  (12)  and 
(14)  in  order  to  obtain  a  more  convenient  expression  of  r\  in 
terms  of  0,  Bim  and  the  surface  concentration  of  the  reference 
species 
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2.3.  Reactor 

Mass,  energy  and  momentum  balances  are 

dC  _  at]Ry(  1  -£b)  C  f  1  dP  P 
d  z  uz 


d  T 


Cj  \RqT  d z  RgT2  dz 
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Fig.  1.  Simulation  results  and  industrial  data.  CO  conversion  (x)  vs.  dimen¬ 
sionless  reactor  length  (z). 


Table  1 

Basic  operating  conditions 

T  (K) 

473 

P  (atm) 

5 

F  (mol  s-1) 

0.32 

Rp  (cm) 

0.285 

V  (cm3) 

3927 

Jco 

0.06 

>’w 

0.40 

yco2 

0.12 

yn2 

0.40 

ycn4 

0.02 

dr  =  (-AH^R^l  -  eb)  A 
dz  Cpm 


(18) 


dr 

dz 


=  9869  x  10“7 
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3.  Algorithms 


Fig.  2.  Dimensionless  reactor  volume  (t)  vs.  CO  outlet  concentration  [ppm] 
at  different  P  [atm]. 


Simple  first  order  differential  equations  are  solved  by  the 
fourth  order  Runge-Kutta  method. 

However,  non-linear  boundary  value  problems  cannot  be 
solved  as  simply  as  first  order  initial  value  problems.  The 
inverse  shooting  method  was  used  to  accomplish  this  task 
[37,38]  (see  Appendix  A). 

When  this  method  is  applied,  the  order  of  integration  is 
reversed;  therefore,  the  origin  of  coordinates  corresponds  to 
the  pellet  surface,  while  the  end  is  the  centre  of  the  pellet. 

3. 1.  Validation  of  the  model 

The  model  was  compared  with  experimental  data  from  two 
industrial  reactors  reported  by  Elnashaie  and  Elshishini  [39]. 
The  kinetics  was  taken  from  Rase  [40] .  Results  are  shown  in 
Fig.  E  As  it  can  be  seen,  the  agreement  is  quite  satisfactory. 

4.  Results  and  discussion 

The  analysis  was  carried  out  considering  that  the  gas  mix¬ 
ture  entering  the  reactor  has  a  composition  similar  to  that  of 
an  ethanol  steam  reformer  outlet,  namely,  CO,  CO2,  H2  and 
CH4.  Also,  the  study  was  carried  out  considering  that  the  H2 
produced  must  feed  a  lOkW  PEM  fuel  cell. 

Unless  otherwise  informed,  the  analysis  was  performed 
at  time  0  -  without  deactivation  -  and  considering  adiabatic 
operation. 

Operating  conditions  are  listed  in  Table  1 . 

A  comment  should  be  made  about  the  pressure  value.  As 
the  PEM  fuel  cell  feed  pressure  is  near  4.5  atm  [44],  it  seems 


to  be  more  convenient  to  increase  the  pressure  before  the 
WGS  stage  rather  than  before  the  fuel  cell.  Fig.  2  indicates 
the  effect  of  pressure  on  the  WGS  performance.  It  can  be  seen 
that  the  CO  conversion  is  favoured  when  pressure  increases. 

4.1.  Effect  of  pellet  radius 

Fig.  3  shows  that  CO  profiles  depart  from  1  regardless  of 
the  pellet  radius,  indicating  that  the  bulk  and  surface  com¬ 
positions  can  be  considered  equal.  In  addition,  these  profiles 
become  steeper  as  the  particle  size  increases.  The  mass  trans¬ 
fer  resistance  in  the  film  is  only  detected  when  the  gas  flow 
in  the  reactor  falls  to  approximately  a  5%  of  the  value  given 
in  Table  1  [41]. 

Temperature  profiles  are  flat  throughout  the  pellet  radius 
(Fig.  4),  regardless  of  its  size.  Nevertheless,  each  profile 
departs  from  a  different  value,  clearly  showing  the  effect  of 


Fig.  3.  Dimensionless  CO  concentration  (f)  vs.  dimensionless  pellet  radius 
(oj)  for  different  Rp . 
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Fig.  4.  Dimensionless  temperature  ( 6 )  vs.  dimensionless  pellet  radius  ( co ) 
for  different  Rp . 


Fig.  7.  Dimensionless  reactor  volume  (t)  required  for  a  given  exit  CO  con¬ 
centration  at  different  1°. 


Fig.  5.  Pressure  drop  (A P)  vs.  dimensionless  reactor  length  (z)  for  different 
*P- 

the  film.  When  the  gas  flow  is  decreased,  this  effect  is  more 
pronounced  on  temperature  profiles  than  on  CO  profiles  [41]. 

The  effect  of  the  pellet  size  on  the  pressure  drop  and  CO 
conversion  is  shown  in  Figs.  5  and  6,  respectively.  It  can 
be  seen  that  A P  increases  and  the  CO  conversion  decreases 
when  the  pellet  radius  decreases.  Thus,  the  size  of  the  pellet 
is  a  trade-off  between  pressure  drop  and  conversion. 

4.2.  Effect  of  the  inlet  temperature 

Fig.  7  indicates  the  reactor  volume  required  to  reach  a 
given  CO  exit  concentration.  As  can  be  seen,  on  one  hand,  if 
the  CO  exit  concentration  is  not  critical,  high  inlet  tempera¬ 
tures  are  more  convenient  due  to  the  thermal  effects  on  the 
kinetics.  On  the  other  hand,  for  higher  CO  removal,  low  inlet 
temperatures  are  more  suitable  due  to  the  thermal  effects  on 
the  equilibrium.  The  strong  dependence  of  the  reactor  vol- 


Fig.  8.  Reaction  rate  ( Rv )  vs.  dimensionless  reactor  length  (z)  for  isothermal 
and  adiabatic  operation. 

ume  on  CO  concentration  values  near  equilibrium  should  be 
outlined. 

4.3.  Adiabatic  versus  isothermal 

The  reaction  rate  and  conversion  profiles  along  the  reactor 
for  these  thermal  regimes  are  shown  in  Figs.  8  and  9,  respec¬ 
tively.  The  reaction  rate  for  the  adiabatic  reactor  shows  a 
maximum  because  of  the  opposite  effects  of  concentration 
and  temperature,  whereas  in  the  isothermal  case,  the  reaction 
rate  always  decreases  (Fig.  8).  In  the  first  part  of  the  reactor, 
where  kinetic  effects  are  preponderant,  the  reaction  rate  for 
the  adiabatic  case  is  larger  than  that  in  the  isothermal  case. 
Thus,  profiles  develop  more  rapidly.  However,  as  the  reaction 
is  exothermic,  higher  temperatures  obtained  in  the  adiabatic 
reactor  yield  a  lower  conversion  (Fig.  9). 


Fig.  6.  CO  conversion  (x)  vs.  dimensionless  reactor  length  (z)  for  different  Fig.  9.  CO  conversion  (x)  vs.  dimensionless  reactor  length  (z)  for  isothermal 

Rp  •  and  adiabatic  operation. 
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CO  [ppm] 


Fig.  10.  Dimensionless  reactor  volume  (t)  required  for  a  given  outlet  CO  Fig.  12.  Dimensionless  reactor  volume  (t)  required  for  a  given  outlet  CO 

composition  [ppm]  for  two  different  catalyst  activities.  composition  [ppm]  for  different  operation  times. 


4.4.  Effect  of  activity 

Fig.  10  shows  the  reactor  volume  required  to  reach  a  given 
exit  concentration  using  the  actual  catalyst  ( k=  100%)  and  a 
catalyst  twice  as  active  as  the  original  one  (k  =  200%). 

It  can  be  seen  that  the  reactor  volume  is  not  as  sensitive 
to  the  catalyst  activity  as  it  would  be  expected.  When  the 
activity  was  increased  by  a  100%,  the  average  conversion 
increased  only  a  35%  at  the  inlet,  where  the  reaction  rate  is 
the  highest.  This  is  due  to  the  fast  approach  to  equilibrium 
inside  the  pellet,  which  causes  the  effectiveness  factor  value 
to  decrease.  Thus,  the  pellet  acts  as  a  “buffer”. 


Fig.  13.  Activity  (a)  and  temperature  ( T)  vs.  dimensionless  reactor  length 
(z)  at  a  given  operation  time. 


4. 5.  Effect  of  reactor  diameter 

For  a  given  pellet  size,  pressure  is  the  most  strongly 
affected  variable  by  the  reactor  diameter.  It  can  be  observed 
from  Fig.  1 1  that,  for  diameters  larger  than  7.5  cm,  the  reac¬ 
tor  can  be  considered  isobaric.  The  pressure  drop  is  relevant 
for  diameters  smaller  than  3.5  cm.  It  can  be  concluded  that 
there  is  a  minimum  diameter  given  by  the  pressure  drop  and  a 
maximum  diameter  given  by  the  plug  flow  condition  (in  Eq. 
(1)). 


Fig.  14.  Effectiveness  factor  (ij)  vs.  dimensionless  reactor  length  (z)  at  dif¬ 
ferent  operation  times. 


4.6.  Effect  of  the  operation  time 

In  Fig.  12,  the  reactor  volume  required  to  reach  a  given  CO 
exit  concentration  for  different  operation  times  is  shown.  Due 
to  the  deactivation,  for  a  given  CO  concentration  the  reactor 


Fig.  1 1 .  Pressure  drop  (A P)  vs.  dimensionless  reactor  length  (z)  for  different 
Dr. 


volume  increases  with  the  operation  time.  A  convergence 
of  curves  is  observed  for  t  — >  oo  because  of  the  existence 
of  a  non-zero  residual  activity  (see  Fig.  12).  The  effect  of 
the  deactivation  on  the  reactor  performance  is  similar  to  that 
observed  on  Fig.  10,  since  in  both  analyses  there  is  a  change 
in  the  activity. 

Since  the  activity  only  depends  on  the  temperature 
(besides  time),  a  temperature  profile  produces  an  activity  pro¬ 
file  along  the  reactor.  This  is  shown  in  Fig.  13. 

As  the  activity  decreases,  the  effectiveness  factor  increases 
(Fig.  14)  because  the  pellet  is  more  uniformly  used.  Once 
more,  the  effectiveness  factor  acts  as  a  “buffer”. 

5.  Conclusions 

In  this  work,  the  simulation  of  the  WGS  reactor  was  per¬ 
formed,  applying  the  one-dimensional  heterogeneous  model. 
In  the  studied  cases,  the  isothermicity  of  the  pellet  and  the 
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existence  of  temperature  gradients  in  the  film  were  demon¬ 
strated  (Fig.  4).  This  gradient  increases  with  increasing  pel¬ 
let  size  and  decreasing  gas  flow.  As  regards  mass  transfer, 
concentration  profiles  are  developed  inside  the  pellet,  and 
become  steeper  with  increasing  pellet  size.  On  the  other  hand, 
film  resistance  can  be  neglected  (Fig.  3). 

Although  under  certain  conditions  isothermal  operation 
would  allow  higher  conversions,  these  improvements  do  not 
seem  significant;  thus,  adiabatic  operation  is  chosen  due  to 
its  simplicity. 

The  COPROX  reactor,  which  should  be  downstream  with 
respect  to  the  WGS  stage,  can  handle  an  inlet  gas  mixture  with 
a  CO  concentration  between  2000  and  20000  ppm  [42,43]. 
According  to  Fig.  7,  the  required  reactor  volume  is  strongly 
sensitive  to  the  CO  outlet  concentration  for  CO  levels  lower 
than  6000  ppm  (depending  on  the  inlet  temperature).  For  a 
CO  concentration  of  3580  ppm,  operating  adiabatically  with 
7°  =  200  °C,  P  =  5  atm,  Rv  -  0.285  cm  (commercial  size),  the 
reactor  volume  for  a  lOkW  PEM  fuel  cell  should  be  of 
approximately  41  (Table  1).  A  reactor  diameter  of  7.5  cm 
and  a  length  of  around  90  cm  satisfy  the  plug  flow  hypothe¬ 
sis  and  produce  a  negligible  A P  (see  Figs.  5  and  11).  These 
values  were  obtained  at  t-  0,  namely  a  -  1 .  If  this  analysis 
is  repeated  for  a  CO  concentration  of  7 160  ppm  -  twice  as 
higher  as  the  original  value  -  the  reactor  volume  is  reduced 
by  74%,  that  is  to  say,  to  around  1 1,  with  a  reactor  diameter 
of  4.8  cm  and  a  length  of  57.6  cm,  dimensions  that  fulfil  in 
Eq.  (1). 

The  reactor  size  can  also  be  reduced  by  working  with  a 
smaller  pellet  size  (due  to  the  rise  of  the  effectiveness  factor 
and  residence  time),  always  checking  the  pressure  drop.  If 
the  pellet  size  is  0.1  cm  the  pressure  drop  is  still  very  low 
(see  Fig.  5).  In  this  case,  for  a  CO  outlet  concentration  of 
7160  ppm,  the  reactor  volume  is  reduced  by  90%,  i.e.  0.41. 
Taking  into  account  the  catalyst  deactivation,  the  reactor  vol¬ 
ume  can  be  reduced  by  84%,  namely  to  0.64 1. 

As  a  car  requires  a  60  kW  PEM  fuel  cell,  an  alternative 
would  be  to  place  six  reactors  of  the  same  size  in  paral¬ 
lel.  A  single  reactor  with  a  volume  the  sum  of  six  reactors 
does  not  seem  to  be  an  appropriate  alternative,  as  depicted  in 
Figs.  6  and  12. 

So,  performing  an  overall  analysis,  and  taking  into  account 
the  operation  variables  and  reactor  and  pellet  geometry,  it  is 
possible  to  achieve  reasonable  volumes  to  the  WGS  reactor. 
These  volumes  will  also  depend  on  the  inlet  conditions  that 
a  COPROX  reactor  can  handle  -  which  are  in  fact  the  WGS 
outlets  -  considering  that,  the  lower  the  WGS  conversion  is, 
the  higher  the  amount  of  oxidized  hydrogen  in  the  COPROX 
stage  will  be. 
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Appendix  A 

A.l.  The  shooting  method 

This  method  consists  of  two  stages:  integration  along  the 
spatial  position  and  correction  of  the  trial  value.  The  trial 
value  is  the  dimensionless  surface  concentration  of  the  refer¬ 
ence  species. 

Once  the  reference  species  profile  is  solved  inside  the  pel¬ 
let,  the  effectiveness  factor  can  be  calculated  by  means  of 
Eq.  (16)  and  used  to  compute  Eqs.  (17)  and  (18).  These  two 
equations,  together  with  Eq.  (19),  are  used  to  calculate  bulk 
conditions  throughout  the  reactor. 


A.  2.  Integration 


Making  Eq.  (7)  dimensionless  and  introducing  the  change 
of  variable  co  =  RpR  Q ,  Eq.  (20)  can  be  derived 


V 


-aRpRy  2  d§ 
DefC%  1  —  co  d  co 
— aRp  R\ 

3  DefCg 


T2  >  0 


(20) 


A.  3.  Correction  of  trial  value 


With  the  trial  value  v,  the  boundary  condition  at  the  centre 
of  the  pellet  should  be  verified.  This  condition  is  expressed 
as  an  error  that  must  be  driven  to  zero.  The  error  function 
defined  in  this  way  is  a  function  of  the  trial  value 

e(v)=^-(l,v).  (21) 

000 

The  mathematical  “trick”  consists  in  considering  v  as 
another  independent  variable  (besides  the  position).  Never¬ 
theless,  the  error  function  so  defined  is  only  a  function  of 
v  because  the  concentration  derivative  is  always  evaluated 
at  co  =  1 .  £  is  equal  to  zero  when  the  first  derivative  of  the 
concentration  at  the  centre  is  zero. 

From  other  point  of  view,  it  can  be  said  that  s(v)  =  0  when 
v  is  the  root  of  the  function.  Therefore,  the  situation  is  treated 
as  a  root  finding  problem.  The  Newton-Raphson  method  is 
implemented  to  achieve  this  goal 


(22) 


However,  the  denominator  in  Eq.  (22)  has  not  been  defined 
yet.  Two  new  variables  are  introduced 


9X 

3o) 


(23) 

(24) 
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Combining  Eqs.  (21),  (23)  and  (24),  it  can  be  seen  that 


d  s 

fO,  v)  =  —(v) 
dv 


(25) 


The  differential  equation  that  allows  \/r( l,v)  to  be  calculated 
is  obtained  by  derivating  Eq.  (20)  with  respect  to  v 


(EX 
d  co2 


=  < 


—aRp  dR\ 
PefCS 
—oiRp  dR\ 


X  + 


1  —  co 


3£>efCg 


xj/  Q  >  0 


£2  =  0 


(26) 


This  differential  equation  is  solved  together  with  Eq.  (20). 
The  value  of  xj/  at  co  =  1  is  used  to  compute  Eqs.  (21)  and  (22) 
and  update  the  trial  value. 
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